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Abstract 

We present the numbers of ice model and eight-vertex model configurations (with 
Boltzmann factors equal to one), I{n) and E(n) respectively, on the two-dimensional 
Sierpinski gasket SG{n) at stage n. For the eight-vertex model, the number of con- 
figurations is E{n) = 2^(^"+-^)'^ and the entropy per site, defined as liui^^oo In E{n)/v 
where v is the number of vertices on SG{n), is exactly equal to In 2. For the ice model, 
the upper and lower bounds for the entropy per site lim^^oo ln/(n)/t> are derived in 
terms of the results at a certain stage. As the difference between these bounds converges 
quickly to zero as the calculated stage increases, the numerical value of the entropy can 
be evaluated with more than a hundred significant figures accurate. The corresponding 
result of ice model on the generalized two-dimensional Sierpinski gasket SGb{n) with 
6 = 3 is also obtained. For the generalized vertex model on SG3{n), the number of con- 
figurations is 2(^^^"+'^)/^ and the entropy per site is equal to |ln2. The general upper 
and lower bounds for the entropy per site for arbitrary b are conjectured. 

1 Introduction 

The ice model was introduced by Pauling to study the residual entropy of water ice [T] , and was 
solved exactly by Lieb on the square lattice [^ |3] . The eight- vertex model is a generalization of 
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the ice-type (six- vertex) models [HIS] and solved by Baxter for the zero- field case PE]. On the 
triangular lattice, there are 20 vertex configurations for the ice rule [S] and the 32- vertex model 
was considered |9]. There is a correspondence between such model and the Ising model, while 
there are other related models, see, for example, [10]. It is of interest to consider the ice model 
and eight- vertex model on self-similar fractal lattices which have scaling invariance rather than 
translational invariance. Fractals are geometric structures of non-integer Hausdorff dimension 
realized by repeated construction of an elementary shape on progressively smaller length scales 
[TT| [12] . A well-known example of fractal is the Sierpinski gasket which has been extensively 
studied in several contexts [T3l[Tl[l5l[T6l[I71[l8l[nil20l[ai[22l|23l|2l[25l[26l[^ 
We shall derive the recursion relations for the numbers of ice model and eight-vertex model 
configurations with Boltzmann factors equal to one on the two-dimensional Sierpinski gasket, 
and determine the entropies. We shall also consider the number of ice model configurations 
on a generalized two-dimensional Sierpinski gasket. 

2 Preliminaries 

We first recall some relevant definitions in this section. A connected graph (without loops) 
G = (y,E) is defined by its vertex (site) and edge (bond) sets V and E [301 EI]- Let 
v{G) = \V\ be the number of vertices and e{G) = \E\ the number of edges in G. The 
degree or coordination number ki of a vertex Vi E V is the number of edges attached to 
it. A /c-regular graph is a graph with the property that each of its vertices has the same 
degree k. Let us consider a 4-regular graph first, and assign an orientation on each edge. For 
the ice model, the ice rule must be satisfied at every vertex. Namely, the number of arrows 
pointing inward at each vertex must be two and the number of arrows pointing outward is also 
two. There are six possible different configurations of arrows at each vertex. For the eight- 
vertex model, the number of arrows pointing inward at each vertex must be an even number. 
There are now eight possible different configurations of arrows at each vertex, including the 
six configurations of the ice model plus sink (all four arrows pointing inward) and source 
(all four arrows pointing outward). For the vertex with degree six, the number of arrows 
pointing inward and the number of arrows pointing outward are both three. There are twenty 
different arrow configurations at such vertex. For the 32-vertex model, the number of arrows 
pointing inward at degree-6 vertex must be an odd number. The thirty-two different arrow 
configurations includes the twenty configurations of the ice model plus only one arrow (six 
possibility) pointing outward and only one arrow (six possibility) pointing inward. In general, 
one can associate an energy to the vertex for each configuration. All such Boltzmann weights 
are set to one throughout this paper. 

Let us denote the total number of ice model configurations on a graph G as /(G) and that 



of eight- vertex model configurations as E{G). The entropy per site for the ice model is given 

by 

biG= hm — tt:;^ , (1) 

where G, when used as a subscript in this manner, implicitly refers to the thermodynamic 
limit. Similarly, the corresponding entropy per site for the eight- vertex model is denoted as 
Se,g- We will see that the limits Si^g and Se,g exist for the Sierpinski gasket considered in 
this paper. 

The construction of the two-dimensional Sierpinski gasket SG{n) at stage n is shown in 
Fig. [H At stage n = 0, it is an equilateral triangle; while stage n + 1 is obtained by the 
juxtaposition of three n-stage structures. The two-dimensional Sierpinski gaskets has fractal 
dimensionality D = In 3/ In 2 [13], and the numbers of edges and vertices are given by 

e{SG{n)) = 3"+i , (2) 

v{SG{n)) = ^[3- + l]. (3) 

Except the 3 outmost vertices which have degree 2, all other vertices of SG{n) have degree 4. 
In the large n limit, SG is 4-regular. 





A 

SG{0) SG{1) SG{2) SG{3) 

Figure 1: The first four stages n = 0, 1, 2, 3 of the two-dimensional Sierpinski gasket SG{n). 

The two-dimensional Sierpinski gasket can be generalized, denoted as SGb{n), by intro- 
ducing the side length b which is an integer larger or equal to two [32j. The generalized 
two-dimensional Sierpinski gasket at stage ra + 1 is constructed with b layers of stage n struc- 
tures. The two-dimensional SG^in) with 6 = 3 at stage n = 1, 2 and 6 = 4 at stage n = 1 are 
illustrated in Fig. |2j The ordinary two-dimensional Sierpinski gasket SG{n) corresponds to 
the 6 = 2 case, where the index b is neglected for simplicity. The Hausdorff dimension for SGi, 



is given by D = In ( ^ ) / In 6 [32] • Notice that SG^ is not /c-regular even in the thermodynamic 
hmit. We shall use simplified notations hin) and Eh{n) for the numbers of ice model and 
generalized vertex model configurations on SGb{n). 




SG-i{l) 5^3(2) SG^il) 



Figure 2: The generalized two-dimensional Sierpinski gasket SGb{n) with 6 = 3 at stage n = 1,2 and b ~ 4 
at stage n — 1. 



3 The number of ice model configurations on SG(n 

Denote the number of ice model configurations on the two-dimensional Sierpinski gasket SG{n) 
as /(^). In this section we derive its entropy per site in detail. As the three outmost vertices 
of the Sierpinski gasket have degree two, they are exempt from the ice rule in the calculation 
of I{n). For the two edges connected to each of these outmost vertices, each of them can be 
directed either inward or outward independently. Let us define the quantities to be used. 

Definition 3.1. Consider the generalized two-dimensional Sierpinski gasket SGb{n) at stage 
n. (i) Define gb{n) as the number of ice model configurations such that one certain edge 
connected to an outmost vertex, say the topmost vertex in Fig. 0, is directed inward and the 
other five edges connected to the outmost vertices are directed outward, (ii) Define pab{n) as 
the number of ice model configurations such that the two edges of a certain outmost vertex, say 
the left one in Fig. are directed outward; the two edges of another certain outmost vertex, 
say the right one in Fig. 0, are directed inward; only one of the edges of the third outmost 
vertex is directed inward and it is on the same side of the Sierpinski gasket as one of the edges 
of the second outmost vertex. (Hi) Define pbb{n) as the numbers of ice model configurations 
such that the two edges of a certain outmost vertex, say the left one in Fig. are directed 



outward; the two edges of another certain outmost vertex, say the right one in Fig. 0, are 
directed inward; only one of the edges of the third outmost vertex is directed outward and it is 
on the same side of the Sierpinski gasket as one of the edges of the second outmost vertex, (iv) 
Define pcb{n) as the numbers of ice model configurations such that all three outmost vertices 
have one edge directed inward and one edge directed outward, while two certain directed-inward 
edges are on the same side of the Sierpinski gasket, say the upper-right side in Fig. (v) 
Define pdb{n) as the number of ice model configurations such that all three outmost vertices 
have one edge directed inward and one edge directed outward, while all three directed-inward 
edges are on the different sides of the Sierpinski gasket with a certain direction, say clockwise 
in Fig. 0. (vi) Define rb{n) as the number of ice model configurations such that one certain 
edge connected to an outmost vertex, say the topmost vertex in Fig. is directed outward and 
the other five edges connected to the outmost vertices are directed inward. 

Since we only consider ordinary Sierpinski gasket in this and next sections, we shall use 
the notations g{n), pa{n), pb{n), pc{n), pd{n), and r(n) for simphcity. They are illustrated in 
Fig. [3|, where only the outmost vertices and the directions of the edges connected to them are 
shown. In principle, the edges of the three outmost vertices have other possible directions as 
illustrated in Fig. HI but they do not appear in our consideration as discuss below. Because 
of rotational and reflection symmetries, g{n), pa{n), pb{n), pc{n) and r{n) have multiplicity 
six, while pd{n) have multiplicity two. It is clear that the initial values at stage zero are 
fl'(O) = pa,(0) = ]9c(0) = r(0) = and pb{fS) = pd{0) = 1. For the purpose of obtaining the 
asymptotic behavior of I{n) in this section, g{n) and r{n) are not needed such that 

I{n) = 6pa{n) + 6pb{n) + 6pc{n) + 2pd{n) (4) 

for non- negative integer n. The reason of missing g{n) and r{n) is due to their zero value at 
stage zero. For example, g{n + l) may contain a term like g{n)pb'^{n) in its recursion relations, 
but all such terms are equal to zero since g{0) = 0. However, g{n) and r{n) are nonzero for 
n > in the next section for the eight- vertex model. Now the four quantities pa{n), pb{n), 
pc{n) and pd{n) satisfy recursion relations. 









g{n) pa{n) pb{n) pc{n) pd{n) r{n) 



Figure 3: Illustration for the configurations g{n), pa(n), pb{n), pc{ri)., pd(n), and r{n). Only the three 
outmost vertices and the directions of the edges connected to them are shown explicitly. 








f(n) ha{n) hh{n) hc{n) hd{n) 








qa{n) qh{n) qc{n) qd{n) s{n) 



Figure 4: Illustration for other possible edge directions which are not used here. Out of the six edges 
connected to the three outmost vertices, f{n) has all of them directed outward; ha{n), hb{n), hc(n), hd(n) 
have four of them directed outward; qa{n), qb(n), qc{n), qd{n) have four of them directed inward; s{n) has all 
of them directed inward. 

Lemma 3.1. For any non-negative integer n, 

pa{n + 1) = [pa{n) + ph{nyf'[pa{n) + 3pc(ri) + pd{n)] , (5) 

pb{n + 1) = [pa{n) + ph{nyf'[pb{n) + ?)pc{n) + pd{n)] , (6) 

pc{n + 1) = pa^{n)ph{n) + pb'^ (n)pa{n) + [3pc{n) + pd{n)]'^ , (7) 

pd{n + 1) = pa^{n) + ph^{n) + [?>pc{n) + pd{n)f . (8) 

Proof The Sierpinski gaskets SG{n + 1) is composed of three SG{n) with three pairs of 
vertices identified. At these identified vertices, the ice rule should be satisfied. As illustrated 
in Fig. \5\ the number pa{n + 1) consists of two cases. For the first case, the top SG{n) is 
always counted by pa{n). The identified vertex of the lower two SG{n)^s has four possible 
configurations such that for each of these SG{n)^s one edge is directed inward and the other 
edge is outward as represented by a big circle in Fig. [51 They are counted by pa{n)[pa{n) + 
p6(n)]^. One may wonder why the other two possible configurations satisfying the ice rule are 
missing. If the two edges of the left SG{n) are directed outward and the two edges of the right 
SG{n) are directed inward, then the left SG{n) is counted by ha{n) and the right SG{n) is 
counted by qa{n). However, such terms have no contributions because ha{n), qa{n) are always 
zero for any n by the same reason discussed earlier. For the second case, the directions of 
the edges of the bottom identified vertex is fixed as shown in Fig. [51 There are four possible 
configurations for the upper-left identified vertex multiplying four possible configurations for 
the upper-right identified vertex, and they are counted by [pa{n) + pb{n)Y[3pc{n) + pd{n)]. 
Eq. ([5]) is verified by combining the results of these two cases. 

The number pb{n + 1) is almost the same as pa{n + 1) except that the two edges of the 
topmost vertex change directions. It follows that the top SG{n) is always counted by pb{n) 
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Figure 5: Illustration for the expression of pa{n + 1). The representation of a big circle at an identified 
vertex corresponds to four possible configurations such that for each SG{n) one edge is directed inward and 
the other edge is outward. 

for the first case, while the results for the second case remain the same, so that Eq. ([6]) is 
verified. 

As illustrated in Fig. |6l the number pc{n + 1) consists of three cases. For the first case, 
the directions of the edges of the three identified vertices are fixed and the number is counted 
by pa'^{n)pb{n). Reversing the directions of the edges of the three identified vertices in the 
first case gives the second case, which is counted by pb'^{n)pa{n). Finally for the third case, 
each of the three identified vertex has four possible configurations such that for each related 
SG{n) one edge is directed inward and the other is outward. This number if counted by 
[3pc{n) +pd{n)]^, and Eq. ([7]) is verified. 
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Figure 6: Illustration for the expression oi pc{n + 1) 



The number pd{n + 1) is almost the same as pc{n + 1) except that the two edges of the 



topmost vertex change directions. It follows that the first case is counted by pa 



n] 



the 



second case is counted by pb^{n), while the results for the third case remain the same, so that 
Eq. ([H]) is verified. D 

The values of pa{n), pb{n), pc{n), pd{n), I{n) for small n can be evaluated recursively by 
Eqs. ([S])-(IHD as listed in Tabled! These numbers grow exponentially, and do not have simple 
integer factorizations. To estimate the value of the entropy defined in Eq. ([1]), we need the 
following lemmas. 



Table 1: The first few values oi pa{n), pb{n), pc{n), pd{n), I{n). 



n 





1 


2 


3 


4 


pa{n) 





1 


54 


7,953,309 


152,890,249,552,106,555,312,694 


ph{n) 


1 


2 


63 


8,076,510 


152,921,906,677,033,336,640,655 


pc{n) 





1 


131 


146,761,217 


202,319,214,683,073,568,675,255,835 


pd{n) 


1 


2 


134 


146,770,694 


202,319,214,926,381,958,377,247,254 


I{n) 


8 


28 


1,756 


1,270,287,604 


1,620,388,590,888,580,168,157,749,612 



Let us use the notations a{n) = pa{n)/pb{n), f3{n) = pd{n)/pc{n) and 7(n) = pb{n)/pc{n) 
for n > 1. 



Lemma 3.2. a{n) G (0, 1) and (5{n) > 1 for any positive integer n. Furthermore, 

< 1 - a(n + 1) < ^^ [l - a{n)] , 
a{n)-f{nf[l - a{n)]^ < ^^^ + i) _ ^ < J^fi'^ " "H]' 



32/3(n)3 + 7(n)3 
47(n)2[4 + 7(n)] 



32 



<7 



^ ^^ ^(n)^[3 + 7H + /3(n)] 
(n + 1) < i^ ^ 



(9) 
(10) 

7(n)3[l + a(n)-i] + 64a(n)-2/3(n) - 'v" ' -^ - Ig ^ > 

for n > 1, such that the sequence a{n) increases to one, P{n) decreases to one, and 7(n) 
decreases to zero as n increases. 

Proof By Eqs. ([5]) and ([6]), we have 

pb{n + 1) — pa{n + 1) = [pb{n) — pa{n)] [pa{n) + pb{n)]'^ , (12) 

and pa{n) < pb{n) is established by mathematical induction hypothesis a{n) G (0, 1) and the 
initial value a(l) = 1/2. By Eqs. ([7]) and ([8]), we have 



pd{n + 1) — pc{n + 1) = [pa{n) + pb{n)\ [pb{n) — pain)]^ > 

so that (3{n) > 1 for all n > 1. 

Next using the fact that pa{n) < pb{n) and pc{n) < pd{n), we have 



(13) 



pajn) + 3pc{n) + pdjn) pb{n) - pajn) i{n) r . 

< 1 — a[n + 1) = 1 , , ^ . ^ — r^— t^^ < : — ^^ = —. — [1 — a[n)\ , (14) 



pb{n) + 3pc{n) + pd{n) 



Apc{n) 



8 



and 

B(n+1)-1= [pajn) + pb{n)][pb{n) - pa{n)]'^ 

pa{n)pb{n)[pa{n) + pb{n)] + [3pc{n) + pd{n)]^ 



such that 

a{n)'~f{n)^[l — a{n)]'^ 2pa{n)[pb{n) — pa{n)]'^ 
32/3(n)3 + 7(n)3 ~ 2p6(n)3 + [4:pd{n)f 



</3{n + l)-l 



2pb{n)[pb{n) — pa{n)]'^ ■y{n)^[l — a^n)]"^ 
^ [Apc{n)f " 32 ■ ^ ' 



Finial, by Eqs. (|6]) and ([7]), we have 



2^^ [a(n) + l]2[7(n) + 3 + /3(n)] 



so that 



7(n + l) = 7(n) X I rq I /Q^ M3 ' *^^'') 

7(n)'^a(n)^[l + a(n) -^J + [3 + p(n)J'^ 



4a(n)27(n)2 [4 + ^{n)] ^{nY [S + ^{n) + /3(n)] 

< 7(n + 1) < — . (18) 



7(n)3a(n)2[l + a(n)-i] +64/3(n)3 '^ ^ 16 

By Eqs. dH]), (HEl) and (UHl), we have 



l-a{n + l) -f{n) 7(n + 1) 7(71) [3 + 7(n) + /3(n)] 

1 — a(n) 4 ' 7(ri) 16 

(3{n + 2)-l 7(n + l)3[l-a(r2 + l)]2 7(n)5[3 + 7(72) + /3(n)]3[32/3(n)=^ + 7(n)3] 

/3(n + l)- 1 ^ 32[/3(n + 1) - 1] ^ 22ia(n) ' 

(19) 

With a(l) = 1/2, /3(1) = 2 and 7(1) = 2, it is easy to see that a{n) increases to one, /3(n) 
decreases to one and 7(n) decreases to zero as n increases. D 

In passing, we notice that 7(r2)/[l — a(n)] — t- (0, 00) and [/3(ra) — l]/7(n)^/^ — )• (0, 00) in 
the infinite n hmit. 

Lemma 3.3. The entropy for the number of ice model configurations on SG{n) is bounded: 

\npc{m) + ^\n2 < Si,sG < ^lripd{m) + ^\n2 , (20) 



where m is a positive integer. 



Proof By Eq. (j7]) and Lemma [3.21 we have 

pc{n) > 64pc%n - 1) > 6A[6Apc^{n - 2)f > ... > pc{mf"'"^ x QA^^"'"'-'^'>/^ (21) 

for any m < n. Similarly by Eq. ([8]) and Lemma 13. 2[ we obtain 

pd{n) < 64pd%n - 1) < 64[64pd%n - 2)f < ... < pd{mf"'"' x 64(=^""'"-^)/2 ^ ^22) 
so that 
3''-'" In pc(m) +3 (3""'" -1) In 2 < In pc(n) < lnpd{n) < 3"""^ lnprf(m)+3(3"-'"-l) ln2 . (23) 
By Eqs. ([3]) and (jlj), we have 

ln/(n) 2 ro I ^ PcH , g PKn) , ^ P^W i , 2 In pc;(ri) 



T;(5G(n)) 3(3" + 1) L pt/(n) pci(n) prf(n)J 3(3" + 1) 

By the definition of the entropy in Eq. ([T]) and Lemma 13. 2[ 

In I(n) 



'/,5G 



7i-s>oo v{SG{n)) 

2 , r pc(n) pb(n) pa(n)-, , 2\npd(n) 

lim — In 2 + 6i--H- + 6hV4 + 6H^ + lim 



n^oo3(3" + l) L pd{n) pd{n) pd{n)^ n-^oo 3(3" + 1) 

2\n pd(n) , , 

= lim , ^ [ . 25 

n^oo 3(3" + 1) ^ ' 

The proof is completed using the inequality ( l23l) . D 

The difference between the upper and lower bounds for Si^sG quickly converges to zero as 
m increases, and we have the following proposition. 

Proposition 3.1. The entropy per site for the number of ice model configurations on the 
two-dimensional Sierpinski gasket SG{n) in the large n limit is Sj^sg = 0.515648810655.... 

The numerical value of Sj^sg can be calculated with more than a hundred significant 
figures accurate when m in Eq. (!20l) is equal to eight. It is too lengthy to be included here 
and is available from the authors on request. As the square lattice (sq) also has degree 4, it 
is interesting to compare our result with Lieb's on the square lattice p,i3|, Si^sq = 1 1^ | = 
0.431523108677.... We see that the entropy per site on the two-dimensional Sierpinski gasket 
is significantly larger. 

The upper bound given in Lemma 13.31 can be improved further as follows. 



10 



Lemma 3.4. For any integers n > m > 1, 

64(^"~'"~^)/^pc(m)^"~'" < pc{n) 



< 64(3"-'"-i)/V(m)=^""'" [^^^^ + (1 + 



32 



(5{m) - lx3-|(3"-'"-l)/2 



Proof By Eq. (jTj), we know 

pc{n + 1) = pcinf a{nf'-f{nf + a{n)'-f{nf + (3 + p{n))' 



for any n > 1. By Lemma [3.2[ we have for any positive integers n> m>l 

64pc(n - \f < pc{n) < 6Apc{n - 1)^[^^^ + (l + "M ) ] 



(26) 



(27) 



(28) 



Using the formula a„ = c*-^" '" ^-"^^a^ ™ if aj+i = ca| for j = m, m + 1, ..., n — 1, where c is 
a constant, then Eq. f l26|) is estabhshed. D 



Lemma 3.5. T/ie logarithm of the number of ice model configurations divided by the number 
of vertices on SG{n) is bounded: 



2 In 2 (3-™ -3"") In 4 2 3-^" 
3" + 1 1 + 3-" 3 1 + 3-" 

ln/(ra) 2 

^ v{SG{n)) ^ 3" + l 

3-m_3-n 

H <^ln4 + ln 



Inpc(m) 



1 + 3-" 



.2 + lln 


"/3(m) 


3 , ,-] 


} 


[l , ^^"^ 


4 J 


4- 


1 


7("^)h 
32 J 



2 3-™ , , , 

+ 3 X ^^3_„ lnpc(m) 



(29) 



for positive integers n > m > 1. 

Proof The number of ice model configurations can be rewritten as 
I{n) = 6 \pa{n) + ph{n) + pc{n)\ + 2pd{n) 

= 2pc{n)\3[a{n) + l]-^{n) + 3 + P{n) 



then by Lemma 13. 2[ 



3 1 r 3 

8pc{n) 1 + -a{n)j{n) < I{n) < 8pc{n) [5{n) + -'~i{n) 



(30) 



(31) 
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It follows that 



\nl(n) 



< 



—, 31n2 + lnpc(n) < ,^^, ,, 

3(3" + l)L ^ ^ '\ v{SG{n)) 

2 ( 3 

— — -Y-|31n2 + \npc{n) + ln[/3(m) + -7(m 



By Lemma [3.41 the upper and lower bounds of In pc(n) are given by 

on—m 1 



■(3 In 4) + 3"-™ In pc(m) < In pc(n) 



on—m 1 

< <!31n4 + ln 



^ ' +(1 + -^ 



L 



on—m 1 

< <i31n4 + 31n 



32 

+ 



4 



4 



In 



^3 | + 3n— inpc(m) 



1 + 



32 



+ 3"-™ In pc(m) 



Combining Eqs. fl52]) and ( 15^ . the proof is completed. D 

Taking the limit ra — > oo in Eq. fl2U]) . the entropy is bounded as 



< Si^sG - 


_3-m 


[ln4+^ 


In pc(m 


) 




< 3-™|ln 


[l + ^ 


(m)-li 

4 J 


4- 


1 


32 J 



(32) 



(33) 



(34) 



Using the inequality ln(l + e) < e for any e G (0, 1), we obtain the improved upper bound 



0<^/, 



SG 



ln4 H — In pc{m) 
o 



< 3 



/3(m) — 1 'y{m)^ 



96 



(35) 



According to Lemma 13.21 it is clear that the upper bound converges to zero quickly as m 
increases. 



4 The number of eight- vertex model configurations on 

SG{n) 



In this section, we derive the number of eight- vertex model configurations, denoted as E{n), 
on SG{n) and the entropy per site. Now directions of the four edges of a vertex (except the 
three outmost vertices) are allowed to be all inward or all outward. The configurations g{n) 
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and r{n) shown in Fig. Elare involved as their recursion relations contain nonzero terms such 
as pb^{n — 1). We have 

E{n) = Qg{n) + 6pa{n) + 6pb{n) + 6pc{n) + 2pd{n) + 6r(n) (36) 

for non-negative integer n. The recursion relations are lengthy and given in the appendix. As 
the initial values at stage zero are g{0) = pa{0) = pc{0) = r(0) = and pb{0) = pd{0) = 1, it 
turns out that their values are the same for any n > 0: 

g{n) = pa{n) = ph{n) = pc{n) = pd{n) = r{n) = 2*-^" ^'^)/2 _ (^^f) 

We have the following proposition. 

Proposition 4.1. The number of eight-vertex model configurations on the two-dimensional 
Sierpinski gasket SG{n) is E{n) = 2^^^"^^'''^ and the entropy per site in the large n limit is 
Se,sg = ln2. 

Compare again with the square lattice (sq). The special case of Baxter's result [7] with 
Boltzmann weights equal to one also gives SE,sq = In 2. Namely, the entropies per site for the 
two-dimensional Sierpinski gasket and the square lattice are the same. 

5 The number of ice model configurations on SG})(n) 
with 6 = 3 

In this section, we consider the generalized two-dimensional Sierpinski gasket SGb{n) with the 
number of layers b equal to three. For SG^^n), the numbers of edges and vertices are given 

by 

eiSG^in)) = 3 X e'^ , (38) 

7 X fi" -I- 8 

v{SG,{n)) = -—- , (39) 

5 

where the three outmost vertices have degree two. There are (6" — l)/5 vertices of SGsln) 
with degree six and 6(6" — l)/5 vertices with degree four. For each of the vertices with degree 
six, we allows three edges directed inward and the other three directed outward, just like the 
consideration in 20- vertex triangular ice-rule problem [S]. 

By Definition 13. H the number of ice model configurations is Isin) = Gpa^ln) + Gpb^ln) -\- 
Gpc^^n) -\- 2pd^{n). The initial values are the same as for SG: pa^{Q) = pcs^O) = and 
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^63(0) = pd^iO) = 1. Again, (73(72) and r^ln) are zero for any nonnegative n. The recursion 
relations are lengthy and given in the appendix. Some values oi pa^i^n), pb^i^n), pc^{n), pd^i^n), 
Iz{n) are listed in Tabled These numbers grow exponentially, and do not have simple integer 
factorizations. 



Table 2: The first few values of pa3(n), pbs{n), pc^{n), pds{n), hin). 



n 





1 


2 


3 


pa^{n) 





11 


29,665,405,536 


8,329,624,787,357,979,293,987,412,541,852,867,738,187,867,580,946,289,512,122,074,041,155,584 


Ph{n) 


1 


15 


29,990,772,448 


8,329,642,677,826,723,066,417,765,699,958,803,959,673,796,982,684,471,740,765,855,746,621,440 


pc3{n) 





15 


527,746,306,872 


708,045,663,245,136,812,838,888,349,048,042,396,698,172,710,195,432,765,544,715,130,291,741,523,968 


pdsin) 


1 


18 


527,789,051,704 


708,045,663,245,233,047,276,563,406,582,556,247,659,516,105,436,007,524,358,474,481,355,710,267,392 


h{n) 


8 


282 


4,579,993,012,544 


5,664,465,261,566,078,079,800,619,338,522,817,745,538,255,642,031,993,426,552,757,072,040,596,340,736 



By a similar argument as Lemma [3. 3 [ the entropy for the number of ice model configura- 
tions on SG^in) is bounded: 



5 



Tnpc3(m) + 
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\n2 <Sj,sGs< 



7x6™ ^ -^ ' 7x6™ ~ ^'"-^-7x6™ 
where ?7i is a positive integer. We have the following proposition. 



5 15 

In pd^lm) + z — In 2 



7x6^ 



(40) 



Proposition 5.1. The entropy per site for the number of ice model configurations on the gen- 
eralized two-dimensional Sierpinski gasket SG^in) in the largen limit is Sj^sGs = 0.576812423363. 

The convergence of the upper and lower bounds remains quick. More than a hundred 
significant figures for Sj^sGs can be obtained when m in Eq. (BOj) is equal to five. 

6 The number of generalized vertex model configura- 
tions on SGb(n) with 6 = 3 

In this section, we derive the number of generalized vertex model configurations, denoted as 
i?3(n), on SG^ln) and the entropy per site. For the vertices with degree four, the number of 
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arrows pointing inward is even just like the eight-vertex model on the square lattice; while 
for the vertices with degree six, the number of arrows pointing inward is odd just like the 
32- vertex model on the triangular lattice. The configurations gsln) and r3{n) shown in Fig. 
|3]are again involved. We have 

Esin) = Ggsin) + 6pa3{n) + Qpbs{n) + 6pcs{n) + 2pd^{n) + 6r3(n) (41) 

for non-negative integer n. The recursion relations are too lengthy to be included here. 
They are available from the authors on request. The initial values are the same as for SG: 
5'3(0) = ^^3(0) = pc3(0) = r3(0) = and ^063(0) = ^^3(0) = 1. It turns out that their values 
are again the same for any n > 0: 

g^{n) = pas{n) = ph{n) = pc,{n) = pd,{n) = r^in) = 2(8x6"-i8)/5 . (42) 

We have the following proposition. 

Proposition 6.1. The number of generalized vertex model configurations on the generalized 
two-dimensional Sierpinski gasket SG^{n) is E^{n) = 2(^x6"+'^)/^ and the entropy per site in 
the large n limit is Se,sg3 — f lii2. 

From Propositions 14.11 and 16. H we conjecture that the entropy per site SE,sGt ^^r the 
generalized vertex model on the generalized two-dimensional Sierpinski gasket SGb{n) is equal 
to c(6) In 2, where the constant c{b) depends on b and is larger than one for b > 2. 

7 Discussion of ice model configurations 

For the generalized two-dimensional Sierpinski gasket SGi,{n), the numbers of edges and ver- 
tices are given by 

e{SG,{n))=3\^-^^^T , (43) 



6 + 4 



6 + 2 



2 

6(6+1)]" 2(6 + 1 



viSG.in)) = — — ^— ^ + ) ' . (44) 



6 + 2 



2 

The bounds of the entropies for the ice model on SG{n) and SG^i^n) given in sections [3] and 
[5] lead to the following conjecture for general SGhiji). 
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Conjecture 7.1. The entropy per site for the number of ice model configurations on the 
generalized two-dimensional Sierpinski gasket SGb is bounded: 

6 + 2 b + 2 
r77-T^[lnpcfe(m) + 31n2l < Si sg, < ,,,^-,, \\npdb(m) + 3\n2\ , (45) 

where m is a positive integer. 

However, to calculate pcb{m) and pdb{m) for general b may be difficult. We notice that 
Si,sG3 is ^ tiit larger than Si^sg- K is expected that the value of Sj^sg^ increases slightly as b 
increases for the generalized two-dimensional Sierpinski gasket. 
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A Recursion relations for the eight-vertex model on 

SG(n] 



We give the recursion relations for the eight-vertex model on the two-dimensional Sierpinski 
gasket SG{n) here. We will use the simplified notation Qn+i to denote g{n + 1) and similar 
notations for other quantities. For any non-negative integer n, we have 

9n+i = ^gl + ^glpan + ^glphn + ISfi-^pCn + Aglpd^ + ^gnPal + ^gnPhl + lOgnPanPbn 
+12gnpanPCn + Ag^pa^pdn + l2gnpbnPCn + Agnpbnpdn + 4:glrn + pa^ + pbl 
+3palpbn + Spalpcn + pa^pdn + Spblpan + SpblpCn + pb^pdn + QpanpbnPCn 
+2panpbnpdn + AgnpanTn + AgnpbnTn + palrn + p6^r„ + 2panpbnrn , (46) 

pan+1 = pbn+l 

= 2glpan + 2glpbn + "ignpa^ + "igupb^ + GgnPdnPbn + Ggnpanpcn + 2gnpanpdn + GgnpbnPCn 

+2gnpbnpdn + ^glvn + pO^n + P^n + ^P^lpK + ^PalpCn + pa^nP^n + ^pblpttn + SpblpCn 
+pblpdn + GpanpbnPCn + 2panpbnpdn + SgnPttnTn + ^gnPbnTn + 125f„pC„r„ + AgnpdnTn 
+3palrn + SpblVn + GpUnPbnrn + QpanPCnTn + 2panpdnrn + QpbnPCnTn + 2pbnpdnrn 

+4gnrl + 2pa„rl + 2pbnrl , (47) 
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PCn+1 = pdn+1 

= 9I + ^glpa-n + ^glpbn + ^9nPal + Sgnpbl + Qgnpanpbn + 'iglr-n + pal + phi + 27^4 
+pdl + Spalpbn + Spblpan + 27pclpdn + 9pdlpCn + QgnPanrn + Qgnpbnrn + 'ipalr-n 
+3pblrn + QpanpbnTn + 35f„r^ + 3j9a„r^ + 3p6„r^ + rl , (48) 

Tn+i = gnpal + gnpbl + 2gnpanpbn + pal + pbl + Spalpbn + Spalpcn + palpdn + Spblpan 

+3pblpCn + pblpdn + GpanpbnPCn + 2panpbnpdn + '^gnPanTn + '^gnPbnTn + 5pa^r„ 
+5pblrn + lOpanpbnTn + l^panPCnVn + ApanpdnTn + 12p6„]9C„r„ + AphnpdnTn 

+Agnrl + Spanvl + Spfo^r^ + I2pcnrl + Apdnvl + 4r3 . (49) 

B Recursion relations for the ice model on SGsin 



We give the recursion relations for the ice model on the generalized two-dimensional Sierpinski 
gasket SGs{n) here. Since the subscript is 6 = 3 for all the quantities throughout this section, 
we will use the simplified notation pa„+i to denote pas{n + 1) and similar notations for other 
quantities. For any non-negative integer n, we have 

pa„+i = 3pal + pbl + IGpalpbn + GpalpCn + 2palpdn + Spblpan + Sbpa^pbl + bApa^pcl 
+Qpalpdl + 25pblpal + ^Apb^pcl + Gpb^pdl + AOpalpbl + lOSpalpcl + Apalpdl 
+ 108pblpcl + ^pblpdl + 2ApalpbnPCn + Spalpbnpdn + 36palpCnpdn + epb^pa^pCn 
+2pblpanpdn + 36pblpcnpdn + 36palpblpcn + 12palpblpdn + 216palpclpbn 
+ 108j9a„pc„p(i„ + 24panpdnpbn + SQpa^pd^pCn + 24pbnpanPCn + Spb^pa^pdn 
+2lQpblpclpan + lOSpblpclpdn + 2Apblpdlpan + SGpblpdlpCn + 32Apclpalpbn 
+324pclpblpan + Updlpalpbn + Updlpblpan + 324palpblpcl + 36palpblpdl 
+144palpbnPCnpdn + lAApblpanPCnpdn + 216palpblpcnpdn + 324palpclpbnpdn 
+ 108palpdlpbnPCn + 324pblpclpanpdn + lOSpblpdlpanpCn , (50) 

pbn+i = pal + 3pbl + Spalpbn + IQpblpan + QpblpCn + 2pblpdn + 25palpbl + hApalpcl 
+6palpdl + 35pblpal + hApblpcl + Qpb^pdl + AOpalpbl + lOSpalpcl + Apalpdl 
+108pblpcl + Apblpdl + QpalpbnPCn + 2palpbnpdn + 36palpCnpdn + 2ApblpanpCn 
+Spblpanpdn + 36pblpcnpdn + 2ApalpblpCn + Spalpblpd^ + 216palpclpbn 
+108palpclpdn + 2Apalpdlpbn + 36palpdlpCn + 36pblpalpCn + 12pblpalpdn 
+216pblpclpan + lOSpblpclpdn + 2Apblpdlpan + 36pblpdlpCn + 32Apclpalpbn 
+32Apclpblpan + Updlpalpbn + Updlpblpan + 32Apalpblpcl + 3Qpalpblpdl 
+lAApalpbnPCnpdn + lAApblpanPCnpdn + 216palpblpcnpdn + 32Apalpclpbnpdn 
+108palpdlpbnPCn + 32Apblpc\^panpdn + lOSpblpdlpa^pCn , (51) 
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pCn+1 = 5832pcl + 8pdl + pa^phn + ^pa^pCn + pa^pd^ + pb^pa^ + Spb^pCn + pb^pdn 
+ 11664j9c^rf„ + lUpdlpCn + 4.pa^nVbl + "^pb^pal + 9720pc^pci2 ^ io80p<p4 
+Qpalpbl + lQ2palpcl + 6pa^prf^ + 162p63pc3 + Qpblpdl + 4320pc^pd^ 
+21pa^p6„pc„ + 7paf^pbnpdn + 21p6^j9a„pc„ + 7pb'^panpdn + 48palpblpCn 
+ lQpa^pbj^pdn + lQ2panPCnpdn + Mpa^pd^pCn + ASpb^pa^pCn + IGpb^pa^pdn 
+ lQ2pblpclpdn + 5Apbf^pdlpCn + ASGpcl^palpbn + ASGpclpb^pan + ISpdl^pa^pbn 
+18pdlpblpan + 486palpclpbnpdn + 162palpdlpbnPCn + A8Qpblpclpanpdn 
+162pblpdlpanpcn , (52) 

pdn+i = pal + pbl + 5832pc^ + Spd^ + 3palpbn + 9palpCn + Spa^pdn + Spb^pan + dpb^pCn 
+3pblpdn + lieeApclpdn + lUpdlpCn + Spa^pbl + Spb^pal + 9720pc^^pdl 
+1080p<p4 + 2palpbl + lQ2palpcl + Gpalpdl + lQ2pblpcl + Qpblpdl 
+4320pclpdl + 27pa^pbnPCn + dpa^pbnpdn + 27pb'^panpCn + dpb^panpdn 
+36palpblpCn + I2palpblpdn + lQ2palpclpdn + bApaf^pdlpCn + SQpblpalpCn 
+I2pb%palpdn + lQ2pblpclpdn + bApbl^pd^pCn + ASQpcl^pa^pbn + 486pc^p6^pa„ 
+18pdlpa1pbn + ISpdl^pb^pan + ASGpalpc^pbnpdn + lQ2palpdlpbnPCn 
+486pblpclpanpdn + lQ2pblpdlpanPCn ■ (53) 
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